Superconductor-insulator transition driven by local dephasing 
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We consider a system where localized bound electron pairs form an array of " Andreev"-like 
scattering centers and are coupled to a fermionic subsystem of uncorrelated electrons. By means 
of a path-integral approach, which describes the bound electron pairs within a coherent pseudospin 
representation, we derive and analyze the effective action for the collective phase modes which 
arise from the coupling between the two subsystems once the fermionic degrees of freedom are 
integrated out. This effective action has features of a quantum phase model in the presence of a 
Berry phase term and exhibits a coupling to a field which describes at the same time the fluctuations 
of density of the bound pairs and those of the amplitude of the fermion pairs. Due to the competition 
between the local and the hopping induced non-local phase dynamics it is possible, by tuning the 
exchange coupling or the density of the bound pairs, to trigger a transition from a phase ordered 
superconducting to a phase disordered insulating state. We discuss the different mechanisms which 
control this occurrence and the eventual destruction of phase coherence both in the weak and strong 
coupling limit. 

PACS numbers: 74.20.-z, 74.25.-q, 74.20.Mn 



I. INTRODUCTION 

The problem of interacting Cooper pairs and/or 
bosons, together with the possibility of a quantum 
control of long range phase coherence in low dimen- 
sional systems has received considerable interest since the 
experimentsi on homogeneous lead and bismuth films, 
which exhibit a transition from a superconducting to an 
insulating phase as a function of thickness. 

It is known, since the work of Abrahams et ali^, that 
no true metallic behavior can be expected to be observed 
for 2D non-interacting electrons at T=0, because all the 
states will be localized by an arbitrarily small amount of 
disorder. When one includes Coulomb interaction, the 
situation is less clear, but the common belief is that a 
metallic phase should still not appear at T=0 in the pres- 
ence of disorder - though no rigorous proof is available^. 
Yet, in the presence of attractive interactions, one ex- 
pects a superconducting state both at T=0 and finite 
temperature, even in presence of a finite amount of dis- 
order, due to its relative ineffectiveness for a transition 
of Kosterlitz-Thouless type (the degree of relevance being 
given by the Harris criterion-). 

The existence of a superconducting state at T=0 in 
2D systems is then considered to be directly linked to 
that of an insulating state, with no intermediate metallic 
phase present. One hence should be able to observe a 
direct superconductor-insulator transition (SIT) at zero 
temperature in 2D as a function of disorder, interaction 
strength, magnetic field or any other external parameter 
which can drive the system away from the superconduct- 
ing phase. 

There are different theoretical scenarios which are com- 
monly discussed in connection with the salient features 
of such quantum phase transitions: i) dissipative mod- 



els, considering a network of Josephson coupled super- 
conducting grains shunted by resistors^iSiLSi24i&, where 
the competition between the Josephson intergrain cou- 
pling and the charging energy yields an increase of the 
phase fluctuations of the superconducting order param- 
eter and hence leads to a phase disordered state. Dis- 
sipation thereby plays the role of suppressing quantum 
phase fluctuations and thus competes with the charging 
mechanism, ii) Bose-Hubbard models, where the super- 
conducting phase is due to the Bose-Einstein condensa- 
tion of charge-2e bosons and the insulating phase due to a 
proliferation of vortices and localization of pairsii*i2iiii4. 

While the SIT was long thought to be a paradigm 
of the above theories, recent experiments have thrown 
some doubt about that. They revealed, what seems to 
be a low-temperature metallic state which is intertwined 
in such a SIT and thus requires a more appropriate 
theoretical description. In magnetic field tuned experi- 
ments on Mo-Ge samplesi^, in granular superconductors 
(Ga films^- and Pb films',), and in Josephson junction 
arrays^^ a metallic phase has been observed. Moreover, 
the recent experiments of Kapitulnik et alti^, in which a 
metallic phase has been observed sandwiched between 
the superconductor and the insulating phase, suggest 
that perhaps two phase transitions accompany the loss 
of phase coherence in 2D superconductors : (i) a super- 
conductor to a "Bose metal" and (ii) a "Bose metal" to 
an insulator. A "Bose metal" in this context, is thought 
of just a gapless non-superfluid liquid with metallic like 
transport. 

These experimental results have led to reconsider the 
whole issue of the SIT. One recent proposal to handle 
such a new viewpoint of the SIT has been to reexam- 
ine the standard on-site charging model. By includ- 
ing nearest neighbor charging terms it was shown that 
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the resulting uniform Bose metal state lacks any trace 
of either phase or charge order—, due to a competi- 
tion between the order parameters which describe the 
onset of charge order and that of phase coherence. A 
different point of view22i is, that in the quantum dis- 
ordered regime a cancellation happens between the ex- 
ponentially long quasi-particle scattering time and the 
exponentially small quasi-particle population, which ul- 
timately leads to a finite dc conductivity. Finally, intrin- 
sic as well as extrinsic sources of dissipation have been 
suggested as potentially relevant for the occurrence of 
a non-superfluid metallic phase in proximity to the su- 
perconducting phaseiSiSiSi. In this context, the dissi- 
pation is a relevant perturbationSi, which, depending on 
the strength of the coupling to the dissipative source, 
can drive the system from a superconductor-to-insulator 
transition to a superconductor-to-metal transitior^^^^. 
Several other ideas supported a scenario where the sys- 
tem could break into superconducting and insulating " is- 
lands" weakly linked via percolating paths before going 
into a metallic phase, dominated by vortex dissipation^S'. 
This " puddle" scenario matches with the view to describe 
the 2D superconductor-metal transition via supercon- 
ducting islands embedded in thin metal films^i. Detailed 
analysis of the puddle- like model considered strongly fluc- 
tuating superconducting grains embedded in a metallic 
matrix^^, predicting a metal-to- superconductor transition 
with a metallic phase just above the transition domi- 
nated by Andreev reflections between the superconduct- 
ing grains. Finally, recent investigations directed the at- 
tention to new phases with Bose metallic features, where 
the dissipation is dynamically self-generated, as in the 
quantum phase glass model, where disorder in the distri- 
bution of the tunnelling amplitudes and quantum fluctu- 
ations destroy phase coherenceSfi.. 

On similar lines as those dealt with in a great variety 
of such different approaches discussed in the literature, 
our aim here is to investigate a system where the break- 
down of superconductivity situates itself in between the 
case of a "bosonic" and a "fermionic" mechanism for su- 
perconductivity suppression. We focus on determining 
the possible ground states for interacting Cooper pairs, 
in close relation to the classical notions of superfluidity 
and localization of bosons and, on a more general basis, 
how the phase coherence can be tuned to a phase disor- 
dered state whose transport properties may be unconven- 
tional. In particular, we consider within the framework 
of a boson- fermion model (BFM), a system where lo- 
calized bound electron pairs (hard-core bosons) form an 
array of " Andreev" -like scattering centers coupled to a 
fermionic subsystem of itinerant uncorrelated electrons 
(quantum pair-exchange). This scenario goes beyond 
that of pure phase models widely discussed in the lit- 
erature, since here one is dealing with bosonic degrees of 
freedom (localized bound pairs) as well as fermionic ones 
(itinerant electrons) for which, due to the emergence of 
pair correlations, one is dealing with both amplitude as 
well as phase modes. The possibility to tune from short 



to long-range phase coherence arises from the following 
competing effects: 

(i) on the one hand, the short range interaction be- 
tween bosons and fermion pairs (holes) induces a local 
phase locking in a configuration with a quantum super- 
position between bosons and electron pairs, leaving the 
common phase undetermined. 

(ii) on the other hand, the itinerancy of the electrons 
tends to lock and rigidly extend these initially arbitrary 
local phases. As a result, phase coherence develops over 
longer distances by suppressing the quantum fluctuations 
of the local phase, which will involve the dynamics of 
amplitude fluctuations. 

Out of this competition one can recover either a super- 
fluid state in the regime of a small scattering rate or a 
phase disordered state in the limit where the pair ex- 
change dominates and the local quantum phase fluctua- 
tions do not allow for long-range phase coherence. As we 
shall show, the physics described by this scenario strongly 
depends on the concentration of fermions and bosons, 
on the coupling strength, as well as on temperature and 
magnetic fields. The purpose of this study is to analyze 
how a transition from a superconductor to an insulator 
may occur and discuss the features which can give rise 
to unconventional dynamics and, eventually, unconven- 
tional transport properties in the proximity of such a 
transition. 

The outline of the paper is as follows. In Sect. II we 
will introduce the BFM and describe its main qualitative 
features and the phenomena which can be described by 
it. In Sect. Ill, we shall develop a path integral represen- 
tation of this model and derive an effective coarse grained 
action, which, after integrating out the fermionic degrees 
of freedom, is able to describe the low energy dynamics 
of the phase and amplitude modes. In Sect. IV we will 
discuss the phase diagram based on the simplest approx- 
imation of such a coarse grained effective action in terms 
of a phase-only action and explore the transition from 
a superconductor to a phase disordered state as a func- 
tion of the coupling and the density of bosons. In Sect. 
V, we will study the physical features which arise from 
the intrinsic Berry phase term present in our effective ac- 
tion and which arises from the hard core nature of the 
bosons, represented by quantum pseudo spin-i variables. 
In section VI we compare the salient features of this BFM 
scenario with similar scenarios, such as the negative-U 
Hubbard model, the Bose-Hubbard model and Joseph- 
son junction arrays. In the Discussion, section VII, we 
review the main results obtained in this paper and indi- 
cate further developments planed for the near future. 



II. THE BOSON-FERMION MODEL 

The boson-fermion model (BFM) in recent years has 
attracted considerable attention as a model capable of 
capturing basic physical properties in many body sys- 
tems with strong interaction, giving rise to the formation 
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of resonant pair states of bosonic nature inside a reser- 
voir of fermions. Such a scenario was initially proposed 
by one of us (JR) as an alternative to the scenario of 
the hypothetical and yet to be experimentally verified 
bipolaronic superconductivity. It was meant to describe 
the intermediary coupling regime between the adiabatic 
and anti-adiabatic limits in polaronic systems, where an 
exchange between localized bipolarons and pairs of un- 
correlated electrons can be assumed to take place (for 
a recent intuitive justification of such a scenario see for 
instance refi^S). 

As it has turned out, this boson-fermion scenario 
has a much wider range of applicability than that for 
which it was initially proposed and seems to apply to 
very different physical situation such as: hole pairing 
in semiconductors'^*', isospin singlet pairing in nuclear 
matten^, d-wave hole and antiferromagnetic triplet pair- 
ing in the positive-U Hubbard modetf^: (and possibly also 
in the t-J model), entangled atoms in squeezed states in 
molecular Bose Einstein condensation in traps^ and su- 
perfluidity in ultracold fermi gases induced by a Feshbach 
resonance ■ . 

The BFM is reminiscent of an anisotropic Kondo lat- 
tice model in terms of a pseudospin-i, but characterizing 
localized electron pairs instead of localized impurity spins 
as in the Kondo analogue. Its Hamiltonian is given by 

H = (i^-;,)^c+c.. + (AB-2M)^(pf + i) 
+ (cicj- + H.c.) + 5 ^ (p+rr + p-r+) 

The pseudo-spin operators [pf , p~ , pf ] denote the local 
bound electron pairs (bosons) and [t^ — cf.^c^^,T~ = 

CiiCii,Tl = 1 — T^T~] the itinerant pairs of uncorrelated 
electrons. [cf^,Cia] stand for the creation and annihila- 
tion operators of the itinerant electrons (fermions) and 
g is the strength of the boson <^ fermionpair exchange 
interaction. The hopping integral for the itinerant elec- 
trons, which is assumed to be different from zero only 
for nearest neighbor sites, is given by t with a band half 
width equal to D = zt, z denoting the coordination num- 
ber of the underlying lattice. The energy level of the 
bound electron pairs is denoted by A^. The number of 
the ensemble of bosons and fermions being conserved, 
ntot = 'T'FT + "^Fi + 2nB, implies a common chemical 
potential /i for both subsystems, ng, np^^_\^ indicate the 
occupation number per site of the hard core-bosons and 
of the electrons with up and down spin states. 

The exchange coupling between the bosons and the 
fermion pairs can be considered as an effective Andreev- 
like scattering leading to local states which are quantum 
superpositions of the form 

|V';oc)j = y" d0i[cos(6l.,/2) cos(0i)p+ + sin(0j)] x 

[cos(0O +T-+sin(0,/2)sin((/.,)]|O). (1) 



Such states evolve gradually out of the system of local- 
ized dephased bosons and essentially uncorrelated free 
fermions, which characterize the high temperature phase 
of this model, when the temperature is decreased be- 
low a certain T* ~ g where resonant pairing (not bound 
pairs!) starts to be induced in the fermionic subsystem. 
These pair states have already features built in which 
are reminiscent of those which characterize Cooper pair- 
ing of fermions as well as superfluidity of bosons. The 
phases of the two coherent states, corresponding to the 
two subsystems, are the same and hence locked together, 
but are averaged over all angles as a consequence of the 
conserved particle number on any given site, ntot = 2. 
Roughly speaking, the ground state of the system is then 
given by a product state Yii \'<Pioc)i with cosOi — 1 and 
which exhibits no phase correlations on any finite length 
scale. Let us next consider the effect of fermion hop- 
ping between adjacent sites. This will give rise to den- 
sity fluctuations on each of those individual sites and thus 
help to stabilize an arbitrary but finite average value of 
the phases {cjii} over a finite length and time scale. In 
this way the localized bosons and fermion pairs acquire 
itinerancji2^*^Si2i which eventually leads to a superfluid 
state in both subsystems^, provided that the effect of 
the local correlations between the bosons and the fermion 
pairs can be sufficiently diminished, but remaining still 
sufhciently strong to guarantee the formation of pairing 
in the fermionic subsystem. Achieving or not this situa- 
tion will depend on the relative importance of the local 
exchange coupling versus the fermion hopping rate (given 
by the ratio g/t) and as well as on the concentration of 
the bosons, as shown by exact diagonalization study^^ on 
finite clusters of this BFM. 

What we shall attempt in the present study is to de- 
scribe this physics in terms of an effective action for the 
phase and amplitude fields of the bosonic fields. In order 
to achieve this we shall put the discussion on a level which 
is more familiar, namely that one of Josephson junction 
arrays and Bose-Hubbard models. For that purpose let 
us briefly sketch the analogy which exists (up to a certain 
point) between the BFM and those systems which have 
been widely discussed in the literature. A physically pos- 
sible realization of this BFM scenario can be imagined in 
form of a network of superconducting grains embedded 
in a metallic environment and where the only mecha- 
nism of interaction between the grains and the fermionic 
background is that of Andreev reflections. Via such a 
mechanism an electron (hole) is reflected on the grain as 
a hole (electron) leaving behind a surplus of two holes 
(electrons) in the fermionic subsystem and of two elec- 
trons (holes) in the grain. If the grains are such that 
they have a large charging energy, the fluctuations of the 
number of pairs on them are energetically unfavorable 
and hence are largely suppressed. We then have a sit- 
uation where the state of the grains switches essentially 
between zero to double occupancy with respect to the av- 
erage occupation, any time an electron (hole) is reflected 
at the interface of the grain. Thus, the quantum dynam- 
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ics of the single grain can be directly represented by a 
pseudospin-i, in order to account for the doublets which 
represent the two possible states of the grains. 

For such a possible experimental setup, the effective 
sites in the BFM have to be considered as defining a reg- 
ular array of grains and having the same periodicity as 
the underlying lattice on which the fermions move with 
a hopping amplitude t. Moreover the size of the grains 
has to be such that it is much smaller than the distance 
between them. A pictorial view of such an experimental 
set up in given in Fig. ^ The analogy between the BFM 
and the array of superconducting grains which scatter 
pairs of fermions in a metallic matrix via Andreev like 
reflections, may ultimately serve as an experimental de- 
vice on which to test and analyze the theoretical issues 
which will be discussed in this paper. 
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FIG. 1: a) Schematic ID representation of the BFM on a lat- 
tice (top and side view). The bosonic and fermionic particles 
move on two different arrays having the same periodicity: the 
fermions are indicated by circles and the bosons by squares 
on the respective arrays, b) The single site configurations for 
the pseudospin and fermionic variables. 



III. PATH INTEGRAL REPRESENTATION: 
DERIVATION OF THE EFFECTIVE ACTION 

A. Generalities 

Let us now construct an effective action which de- 
scribes the BFM, with the aim to extract the dy- 
namical properties of the low energy degrees of free- 
dom of the phase and amplitude modes for the bosons 
and fermion-pairs. We start by expressing the parti- 
tion function in terms of a coherent-state path inte- 
gral representation^, where the fermionic part is formu- 
lated by means of the usual Grassmann variables and the 
bosonic part is described by a pseudospin-coherent state 



representation^Mi: 



A[*i,*i,6li,0i] 



(2) 



where 



A[i',,^i,6,,(l)^]^ / dTY^[is{l - cos 0i) dr<l>i + 

i 

{Ab - 2Ai)cos(?.] +^*,(r)a-i*.(r). (3) 

T denotes the imaginary Matsubara time variable and 
a Nambu spinor representation for the Grassmann vari- 
ables, related to that of the original fermionic operators 
by 



(4) 



The pseudospin is described by a bosonic field 
which in spherical coordinates is given by Si = 
s{sm 6i cos 4'i, sin 9i sin ipi, cos 6i) (see Fig. 9i de- 

scribes the polar angle of the vector with respect to 
the north pole of the z axis, while 4>i is the azimuthal an- 
gle which defines the angular position of the basal plane 
projection of this vector. The first term of the action A is 
the Wess-Zumino term^i^, ensuring the correct quanti- 
zation of the quantum pseudospin variable. For any path, 
parameterized by 0(t) and 0(t), the contribution of this 
term is equal to i s times the surface area of the sphere be- 
tween this path and the north pole. For closed paths this 
has exactly the form of the Berry phased. The second 
term is linked to the density of the bosons n b (t) through 
the cos 6'(t) dependence of the pseudospin. Finally, the 
last contribution of the action contains the coupling be- 
tween the fermionic and bosonic subsystem through the 
Green's function Gi, , determined by: 



(l* k) G.,{r-r') = 6{T-T') 



(5) 



where Ki = {—dr 

and L = g sin 6'i(T)e* '^•'^'^•'iSy and L* being the conjugate 
field of L. Integrating out the fermionic part, one obtains 
the action in terms of exclusively the bosonic fields: 

A= -TrlnG-^+ ( dr ^[is(I - cos 6',(r)) 9^(/)^(t) -f 

(As -2Ai)cos0,(r)] (6) 

where the trace has to be carried out over all internal as 
well as space-time indices. 

Up to this point no approximation has been made in 
the derivation of the action which describes the coupling 
between the bosonic and fermionic degrees of freedom. It 
is important to notice that the variation of the bosonic 
variable 0{t) describes both, the density fluctuations of 
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6(p(l-cos6) 




FIG. 2: Spherical representation of the pseudospin s including 
the Berry phase factor for one possible trajectory F. The 



Berry phase term is exp[is ?>(1 



where A 



is the area of the surface enclosed in the trajectory F. The 
black part indicate the differential portion of the surface on 
the sphere taken respect to the north pole. 



mr)] 




FIG. 3: Representation of a possible path for the local 
pseudospin motion. The field ^[^(t)] indicates the undula- 
tion of the pseudospin vector along the polar direction as 
it arises from the fluctuations of the average boson density 
{{cos 6)) and of the pairing amplitude {{sin 6)), while it pre- 
cesses around the 2:-axis due to the time evolution of the phase 
variable 0(t)). 



the bosonic subsystem (via the projection of the pseu- 
dospin vector on the z axis, i.e. cos 6'(r), being the lon- 
gitudinal component) and the amplitude fluctuations of 
the fermionic pair field (via the transverse part as pro- 
jection of the pseudospin vector onto the basal plane, i.e. 
sm6{T)) (see Fig. The variable 0(r) determines the 
rotational degrees of freedom of the pseudospin vector, 
expressing its phase dynamics. 

For extracting the relevant terms which control the 
low energy dynamics of the coupled phase and amplitude 
modes and performing an expansion, which is meaning- 
ful in terms of the phase variable, it is judicious to make 
the following steps: i) gauge away the phase dependence 
from the term L which permits to separate the trace 
into a part which does not depend on the phase of the 
bosonic field and another part that contains only spa- 
tial and time variations of (t)^, ii) rewrite the term L 
after the gauge transformation, as a sum of two pieces, 
one not dependent on time (which is linked to the aver- 



age density of bosons) and another term containing the 
fluctuations with respect to its mean value. 

The first operation is performed by applying to the 
operators under the trace the rotation Ui — e"^'('^)°'3/2^ 
where cr,- denote the Pauli matrices. Hence, 



Tr In G- 



(7) 



where 



-dr(7o + i^drcj^iir) - /x]cr3 -I- g sin 6'i(r) CTi 

+ ..p''[</'=(^)-0i(^)]'^3/2 



+ 



(8) 



Since the bosonic density is fixed in average, we now 
separate the part which depends on the polar angle in 
a time independent contribution and its time dependent 
correction. That is, the term sin0i(r) is decomposed into 
its average value (sin^^) (which is determined by fixing 
the density of bosons due to the spherical constraint) plus 
a time dependent contribution 1[9{t)\ which contains the 
fluctuation around its average value (Fig. Thus this 
local field, due to the constraint, will describe both: the 
time dependent variation of the density as well as of the 
pairing amplitude. 

Next, let us write the Green's function in the usual 
form as 



Gj/ - Go,;/ + Sij 

with Go~^ = [—drOQ — ^(73 + g ai] Sij and 
Di + Ki, where 

1 d(j){T) 

2 OT 

^ g l[e{T)]ai 



(9) 



T 



9 = .9(sin[6']). 



(10) 

(11) 
(12) 
(13) 



From this point onward we shall assume the average 
density of bosons to be homogeneously distributed and 
thus given byns = ^(1 + (cos 6*)). This implies that the 
bare coupling g is renormalized to 5 as a consequence 
of the spherical constraint of the pseudospin variable. 
Moreover, since (sind) = \/l ~ {cos 9)^ and the density 
of bosons is fixed in average via a suitable choice of the 
chemical potential and of the on-site bosonic energy, one 
can treat the exchange coupling as an external boson- 
density tunable parameter. The variation in the 9 vari- 
able is then simply related to the variation of the bosonic 
density such that if 9 varies in the range [0, tt] then ns 
varies in the interval [0, 1]. 

Before expanding the trace, let us write down explicitly 
the expression of the zero order Green's function, as it 
will be frequently used in the following steps: 



Go 



Q/3 



Q/3 



{uJn) exp[-i LU„ 



(14) 
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with 



or in a compact form 



,,,2 I ,,,2 



— tUJn —fJ. 



(15) 



and where we introduced ujq — g. 



B. Second order loop expansion 

We now evaluate the contribution of the self energy 
to the effective action. This is done in the usual 
way by making a loop expansion in the traced. We shall 
construct that expansion up to second order in the time 
and space derivatives of the phase variable and in the 
terms which contain both, the fluctuations of the density 
and the amplitude. For that purpose we use the standard 
identity: 

TrlnG^i = Tr In [Gq^^ + S] 



Tr In Gn 



rrln[l + GoS] (16) 



and then expand the second term of this expression up to 
second order in S, such as to keep all the contributions 
up to quadratic order in the gradient of the phase. This 
gives 



Tr In G"^ ^ TrlnGn 



-rr[GoS]-irr[GoS]^ (17) 



The first term of this expression is just a constant and 
does not contribute to the dynamics. In the second term 

Tr[GoS] = Tr[GoT] + Tr[GQD] + Tr[GoK] (18) 

the only parts different from zero are Tr[GoiKi] and 
Tr[GoiDi\. Tr[GoiTij] gives no contribution once one 
makes the trace over the site indices. Tr[GoiDi] intro- 
duces a contribution which is proportional to the chem- 
ical potential multiplied by the time derivative of the 
phase ~ indrcl). We will see below that this contribution 
describes an effective off-set charge (in terms of the ter- 
minology of a similar Josephson junction array scenario) 
and which arises from the total fermionic and bosonic 
static density distribution via their dependence on the 
chemical potential. Tr[GQiKi] describes the lowest order 
fluctuations of the bosonic (fermion pair) density due to 
the presence of the spontaneous pair/hole creation out of 
the condensate. Its direct evaluation gives a contribution 



Al = j dTldT2Trps[G()i{Tl,T2)Ki{T2,Tl 



)] (19) 



where Trps represents the trace over the internal pseu- 
dospin index. With Ki{T2,Ti) = Ki{Ti)S{Ti — T2) and 
integrating over the Matsubara times gives: 



dTTrj,s[Go^{0)K,{T)] 

25Go,^^(0) j dTl[e{T)] 
gtanh[/?LJo] 



Wo 



dTl[6{T)l 



(20) 



g tanh[/?Li;o] 



UJQ 



(21) 



Let us next come to the evaluation of the terms which 
contribute to the quadratic order in this loop expansion 
of the trace. The parts which are non zero in these terms 
are the following: 



A2 = Tr [Go, Go, D,] 

A3 = Tr [Go, Go, K,] 

Ai = Tr [Go, Go, D,] 

As - Tr [Go, T„ Go, T,,] 



(22) 



The term A2 contains processes where the fermionic 
background locally couples at different times to the fluc- 
tuations of the bosonic density and to the phase fluctua- 
tions: 

A2^Tr [Go, Gq, A] = 

J dTidT2Trps[Goi{Ti - T2)Ki{T2)Goi{T2 - ti)A(ti)]- 

(23) 

Since Goi(Ti — T2) depends exclusively on the time dif- 
ferences, we introduce the new variables r = '^^^^ ^-nd 
rj = -^-^^5-^, after which the integral over the Matsubara 
time variables becomes: 



A2 = -y drdr,Trp,[Go,(r)if,(r;-r)Go,(-T)A(r + ?7)]. 

As we are interested in the gradient expansion in the 
bosonic phase and density, related to (/)(t) and 0{t), 
we keep only the lowest order in their time derivatives. 
This is done by considering the expansion of the product 
Ki{r] — t) Diir — rj) around t = in order to separate 
the parts which depend exclusively on the local variable 
of the Green function and those which are linked to the 
phase and the density fluctuations. The expansion then 
reads as follows, 



K^iv) A(ry) + r{ ^A(^) 



dr) 



drj 
0{r') 



(25) 



Due to the symmetry of the Green's functions (even for 
T — r), the linear contribution of the series expansion 
oi Kiij] — t) Di{T — r]) cancels in the effective action after 
integrating over the time. Hence, to lowest order one 
obtains, 

A2 ^ \ l dTdi^Trps[Go^iT)K,ir,)Go,{-T)D,{Ti)] 
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= f / dT[Gol\r)Gol\-r) + Gol\r)Gor{- 
- Gof(T)Go.''(-T)-Gor(T)Gof(-r)] x 



or in a compact form, 



(26) 



(27) 



where the coefficient £^2, after integrating over the Mat- 
subara time, is given by: 



Eo = 



IJ,{(3 (jjQ- sinh [/? cjq]) 
4 wg(l + cosh[/3 Wo]) 



(28) 



Let us next consider the term A^, which expresses 
the coupling of the fermionic field to the fiuctuations of 
bosonic density at different times. In order to extract the 
lowest order gradient contributions, we follow the same 
procedure as that just used above, giving us: 

^3 = 1 J dr dvTrps[Go,iT)Ki{r]-T)Go^{-T)KiiT + r,)] 

(29) 

Performing again the expansion in time up to quadratic 
order in the time derivatives of Kiij), we have: 

in - t) (r, + r) ^ {r,f - . (30) 

and hence can rewrite A3 in the following way: 

A3 = ^JdT dvTrp,[Goi{T)Ki{r])Goi{-T)Ki{r])] + 
\j drdii rrp4Goi(r)K;(7?)Goi(-r)<(77)]. (31) 

Carrying out the trace over the internal indices gives: 
Trj,s[Goi{T)KMGoi{-r)Ki{r))\ = [Gof (r)Gol^(-r) 



+2 Go\\T)Go\\-T) + Go\\T)Gor{-T)] {fl[6{r^)f) 
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(32) 



so that 



A3 = E^^ J drj l[e{r,)]^ - Es, J dr,^^^^\ (33) 

By evaluating the integrals over the Matsubara times we 
obtain: 



E: 



3a 



5^ /3ujog'^ + fi'^ sinh[/3wo] 



2 tjg(l + cosh[/3 Wo]) 



E^h = — 



+ 



gec/i[(/3a;o/2)]^ (-6/3M'a;o + ^3''{-^iM + ^l) 

2 ^ 

6^^ sinh[/3a;o]). 



(34) 



Applying the same procedure for the evaluation of A4 
(which has the same functional form as ^3 but involv- 
ing the coupling between the fermionic degrees of free- 
dom and the phase velocity at different time) one gets 
an analogous expression, providing we discard terms of 
higher order in the derivatives than [^^7-^]^: 



Ai ^IJdT dr,Trps[Goi{T)Di{r,)Go^ - t)A(j?)] 



1„ /• . .5(/.(r?),2 



Co I d,i^r. 



(35) 



Here 



Go = I dT{-[Goj\T)GojH-r) - 2 Gor(T)Gof (-r) + 

Gof (T)Gof (-r)]} , (36) 

which after evaluating the integration over the Matsub- 
ara times gives, 

sech[{Puo/2)]^{Pl^^u;o + g^ sinh[/?a;o]) 
= 2^1 • 

Finally, we come to the evaluation of the last term 
A5 which involves the coupling between phase fluctua- 
tions on different sites and the process of single particle 
hopping. This contribution will yield terms which are 
quadratic in the time derivative of the phase {charging 
like, in the terminology of a similar Josepson junction 
array scenario) and moreover will generate an effective 
hopping induced inter-site phase coupling which turns 
out to be similar to the Josepshon coupling in arrays of 
superconducting grains. 

Considering again the lowest order gradient expansion 
contributions we have: 



A5 — Tr [Goi Tij Goj Tji] 

-IJdT dr]Trp,{Goi{T)Tij{ri)Goji-T)Tji{r])- 



Goiir 



Goj(-t) 



dr] 



(38) 



= Pi + P2 
The terms Pi, 2 are given by: 

Pi =-4/ [^?of(r)Gof(-r)] ^ dr]co^\4>M ' <^M\ 
P2 =-f / drr^ ^^G^^^(r)G^^.i-T)+Gf,(T)Gf,{-T)\ 



dj][- 



,{r]) d(l)j{r))^ 



dr] 



dr] 



(39) 



By carrying out the integration over the Matsubara times 
one obtains for the hopping induced inter-site amplitude 
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and phase coupling: 

Pi ^ Ej J dr]COs[<j),{T])- (j)j{T])] 
f3- (—BbJnsp.nhlBuji 

E.T - 



can be absorbed into the chemical potential and another 
one which is quadratic in Z[6'(t)]. With this, iS^-e can be 
rewritten as: 



tfj 5^ (-/3wosec/i[/?wo/2]^ + 2tanh[/3tJo/2] 



(40) 



Similarly, by calculating the coefficient of the term P2 , 
one obtains the strength of a mutual capacitance term 
(in the terminology of a similar Josephson junction array 
scenario) for neighboring effective sites: 



P2 



drj 



drj 



0/2 



13/2 

exp[/3wo] (-6/3.9^^^0 + P^^litJ-'^ + 



12(l + exp[/3c^o])2u;g 



sinh[/3wo]) 



12(l + exp[^wo])2c^g 



(41) 



The final effective action is then given by the sum of 
all the terms evaluated above which are grouped together 
in form of three different contributions: 



S 



j dr [S^ + S^-e + Sg 



(42) 



S'0, Se, S^-e are the contributions arising from exclu- 
sively (i) the phase dynamics, (ii) the fluctuations of the 
bosonic density and of the amplitude l[ff\ and (iii) the 
coupling between them. They are given by: 



St = Ejia + 'Ci) ( 



(r) 2 



E 

(ij) 



dr 

1^ dMr) dcj^.ir) 



dr 



+ Ej cos[0,;(t) - <Pj{t)] 



+ i 



^0(t) 



E 

i 

E 

i 

i E2 



El l[0^{T)] + E3jmT)f ^ Esbi 



d(l>i{T) 

I s — (1 - cos[6'i(r)]) 



dr ^ 



dr 



mr)] 



(43) 



We expand the Berry phase contribution (the first term 
in S^-o) up to second order in 1[9{t)] and subsequently 
redefine this field as /[^(t)] — a + bl[9{T)], with the time 
independent constants a,b chosen in such a way as to 
eliminate any terms linear in / in the the action S. We 
then find two contributions to the action which are linear 
in d(t>(T) /dr: One which is time independent and which 



E^ 



dr 



(44) 



where ^^(t) is quadratic in Z[0(r)]. 

The effective action thus constructed for the BFM 
(Eqs. 142I43|) ') is similar to that of Josephson junction 
arrays with nearest neighbor, as well as, self capacitance, 
Josephson coupling and off-set charge terms. The action 
for the BFM goes however beyond that for such Joseph- 
son junction arrays in the following respects. We have 
extra terms which control the dynamics of the amplitude 
modes, given by Se and an intrinsic Berry phase term 
which gives rise to a direct phase - amplitude coupling 
(Eq. I44|) . where the dynamical amplitude fluctuations 
would correspond to a time dependent offset charge term 
in an analogous Josephson junction array picture. Fi- 
nally, the Berry phase term, being an intrinsic topologi- 
cal term, may give rise to a Magnus force on a vortex, as 
it will be discussed in section V. 



IV. THE SUPERCONDUCTOR - INSULATOR 
PHASE BOUNDARY 

In a first approach to analyze the stability region of 
long-range phase superconducting coherence, we exam- 
ine the effective action (eqs. I42I43|I . restricting ourselves 
to the phase only part of it. This means a study of the 
SIT driven by a competition between the phase coherence 
induced by pair hopping and the disrupting effect of local 
boson density (or equivalently pair field amplitude) fiuc- 
tuations in the presence of a source term for the bosons 
which controls the global boson density via an effective 
chemical potential. Within such an approximation our 
study is equivalent to that of Josephson junction arrays, 
except that the effective coupling constants entering in 
such an action depend in a highly non trivial way on the 
parameters which characterize the original BFM Hamil- 
tonian. This, as we shall see, will lead to novel features 
concerning the phase diagram with a SIT for the BFM 
when we examine it in terms of the boson-fermion ex- 
change coupling g/t and the boson concentration ns- 

We shall determine the phase diagram by means of the 
so-called coarse graining approximation which has been 
successfully applied for this kind of problem and which 
permits to capture the relevant qualitative and quanti- 
tative features of such a Josephson junction array like 
action^^iiS. It is known that, as a consequence of the 
uncertainty relation between the phase (f)i and the pair 
number operator Qi = i-^, the system can switch from 
a phase ordered to a disordered state. An essential part 
of this study will concern how the relevant parameters 
of the BFM Hamiltonian influence the equivalent am- 
plitudes of the Josephson coupling, the capacitance and 
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the off-set charge terms. Thus, at zero temperature, by 
fixing the bosonic density distribution, we find that the 
variation of the ratio between the Josephson coupling 
energy Ej and the charging energy Ec = which 
controls the phase-density interplay, increases from zero, 
goes through a maximum and then decreases to zero with 
increasing g/t . Or else, if one fixes the coupling g/t, then 
by varying the density of the bosonic distribution one 
is able to control the effective coupling which appears 
in Ej and C.^ bt varying the average bosonic density 
g — 2g^nB{nB — !))■ It is thus immediately evident 
that in the BFM scenario there is a non-trivial interplay 
between the renormalization of the Josephson coupling 
and the charging effect. If one goes to the limit of empty 
{riB = 0) and full bosonic occupation {ub = 1), the 
effective coupling ^ — > 0, and the critical temperature 
consequently reduces to zero. 

The general form of the action wc then have to examine 
is given by: 



5. 



phase 



El I 
9" 



{r) 



- 2 dr - 

COs[(f>.i{T) - <f)j{T)] 



y- 



8 9t 



dr 



(45) 



The first term describes the effect of a static off-set charge 
(qi) The second term contains the physics of the pair hop- 
ping processes, an analog to the Josephson tunnelling 
processes, and has a coupling strength Ej with aij — 1 
if are nearest neighbors and zero otherwise. Fi- 

nally, the third term describes the charging term aris- 
ing from the local exchange of boson and fermion pairs 
and from quasi-particle hopping between nearest neigh- 
bor sites. The strength of this charging type interaction 
is given by Cij = {Co-\-zCi)6ij - Ci 5i,i+p which rep- 
resents an effective general capacitance matrix, with the 
vector p running over the nearest neighbors. Co denotes 
the self- capacitance and Ci the mutual one [z being the 
coordination number). 

As mentioned before, to extract the phase diagram we 
make use of the coarse graining approximation^^. The 
main idea of this approach is to introduce a Hubbard- 
Stratonovich auxiliary field which is conjugate to the av- 
erage of (e"^*) and which plays the role of an order pa- 
rameter for the transition from a superconducting to an 
insulating state. Since the phase transition has a contin- 
uous character, one can expand the action in powers of 
the auxiliary field and determine the occurrence of phase 
coherence by looking at the coefficients of the quadratic- 
term in the limit of long-wavelengths and zero frequency. 

We briefly sketch the main steps of such an approxi- 
mation and adapt it to the present scenario of the BFM. 
The partition function for Sphase is given by the sum of 
all the possible paths of the phase variables in the imag- 



inary time and in the real space: 



Z= /n^'^.exp / Y.-^ 



i d(t>i{T) 



2 dr 

Ej ^ aij cos[0i(r) - 0j(r)] 



dr 



To perform the Hubbard-Stratonovich transforma- 
tion, one rewrites the Josephson coupling term as 
^^ exp[z(/>i]aij exp[j(/)j] and then, by using the usual 
Gaussian identity, introduces an auxiliary field tpi{T). 
The partition function Z then becomes: 



Z = 1 J] DrrD^,Dcl>, exp J j^^ dT[{^) ^ 



^3 



exp 




E 



2 dr 



i(f>i 



ip*e' 



Els- 



dr 



(46) 



Hence, starting from the effective action for the aux- 
iliary field, one can perform an expansion up to second 
order in ip and thus derive the usual Ginzburg Landau 
type free energy functional which permits to determine 
the boundary line between the superconducting and the 
insulating state. The corresponding effective action for 
these auxiliary fields "01 V'* 



Sji, — 



In 



J2{^,{r)e'Mr)_^*(^^y-^Mr)^ 

i 

E- 



dr dr 



(47) 



is then expanded to second order, giving: 



S.. 
where 



^ = , drdr (r, r )^P* (r)V', (r) + 0{^P^), (48) 




,(r,r') = (e4^-«-*^W'])o 



(49) 



denotes the two-time phase correlator, which is equal to 
the second derivative with respect to the auxiliary field 
ip and its conjugate at different time and space positions, 
evaluated around their zero values. By performing the 
functional derivation, one gets the following expression 
for it: 
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x.(./)^ ^n-^t'nr"*''iVr'"'°' <»' 

j ni-D(/)^exp[-S'o] 

with ^0 being the part of the action which contains ex- 
clusively the charging contributions, i.e. : 



So 



exp 



Jo 



1 d(j)i{T] 

2 dr 



(51) 



I g^z(T) d(j)j{T) 

^8 9r dr 



According to the scheme outlined above, we now develop 
the partition function Z up to second order in the auxil- 
iary fields, thus putting it into a familiar form: 



with 



= / drdr ^ipU'^) "^'^(t - t ) - 



(52) 

Mr')- 
(53) 

The determination of the conditions for the boundary line 
between the superconducting and the insulating phase 
then reduces to the explicit evaluation of Xij {t,t ) . The 
result for that has been first obtained in Ref. [48], and 
we sketch below the main steps of this derivation. 

In determining Xij (r, t ) it has turned out to be es- 
sential to treat the phase variable 0i in a compact form. 
In order to separate the imaginary time evolution of the 
phase into a periodic part (jjiir) and an non-periodic part, 
one introduces the following parameterization: 



2'KniT 



(54) 



with (0) — (pi {(3) and rii being an integer which counts 
how many times the phase winds over an angle which is 
a multiple of 27r. 

The use of such a relation allows to express the sum 
over all 0i as an integration over 4>i plus a sum over all the 
possible integer values of the winding number rii. After 
performing a number of suitable algebraic operations^, 
one ends up with the following expression for the two- 
time phase correlator: 



ex 




S!{n,}exp 





^exp 



,(55) 



where Ni = qijl -\- rii. This two-time phase correlator 
is local in space and its time dependence follows an ex- 
ponential behavior, if one assumes that the static off- 
set charge Qi has a distribution which is homogeneous 



in space. In the general case, of a dynamic offset charge 
Qi (r) , this time dependence will be modified and can lead 
to qualitative changes in the nature of the SIT. 

Having obtained the expression for the local two-time 
phase correlator in the time and space representation, 
we can now express the effective Ginzburg-Landau free 
energy functional in the Fourier space in the following 
way: 



F,i, 



1 



n.k 



^ i-'*k{^n) [^k ^ - Xk{^7i)]M^n)\ ■ (56) 



Expanding the inverse matrix Q!,„^ in the form a ^ 



(l/z) + P(aVz2) + . 
is written as: 



this free-energy functional finally 



-y 



2 



Xo ■ 



|V'fe('^n)l'' 



(57) 



The transition line is then given by the condition that 
the coefficient of the quadratic term vanishes in the limit 
of vanishing k and cj, that is: 



1 



zEi 



-Xo(0) = 0. 



(58) 



For a quantitative analysis, one has to determine the 
explicit zero frequency limit of the two-time phase cor- 
relator. As mentioned above, the inclusion of the time 
dependent offset-charge coming from the fluctuations of 
the bosonic density can modify the low frequency behav- 
ior of the local phase correlations. 

In order to get a first insight into the underlying 
physics at play here we evaluate the two time phase cor- 
relator in the limit of a purely local capacitance (the so 
called self-charging limit), by keeping only the on-site 
part in the original structure of the capacitance matrix. 
Under those conditions the evaluation of the two time 
phase correlator at finite frequency gives: 



3^ En 



1 



(59) 



with F[n,] = exp[- Y., WC^^nj] and Zo = E{„.} F[n^]■ 
With these expressions we finally can cast Eq. [SSI in 
the form 



zEj 



exp[-4/3gcn 
l-in'-' 



4i^cE™exph4i;c/3m2] 



0, 



(60) 



expressed in terms of the Josephson coupling Ej and 
the charging energy Ec = G^^/2, as given by the 

Eqs.JSZHOEIJ. 

Given this defining equation for the boundary between 
the superconducting and the insulating phase, we want 
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to see now how the intrinsic dependence of the Joseph- 
son and charging energy on the exchange couphng and 
the density of bosons, manifests itself in the competition 
between phase and boson density degrees of freedom. In 
the loop expansion given in the previous section for the 
BFM, we have obtained the amplitude of the intersite 
phase coupling and the charging effect as a function of 
temperature, the ratio g/t and the bosonic density. It 
is the non trivial dependence of those coupling constants 
on the parameters of the original microscopic Hamilto- 
nian of the BFM which leads to the interesting features 
in the competition between the phase and boson density 
fluctuations. Thus, upon increasing the pair exchange 
coupling it is expected that the Josephson coupling is re- 
duced while the charging energy is increased due to the 
local transfer between fermions pairs and bosons as well 
as due to the single particle hopping of fermions between 
nearest neighbor sites. 

We determine the critical line for two different cases in 
order to highlight the role played by, on the one hand, 
the coupling and, on the other hand, by the effective 
boson density. In Fig. 0] we illustrate the transition 
line T^, separating a phase coherent state from a phase 
disordered one, as a function of the coupling strength 
and the effective boson density ub- The evolution of 
the transition line is non monotonic as a function of g/t 
and goes through a maximum at gmax ~ t. The quantum 
critical point, where the SIT occurs, is given by gcrit ^ 2i. 
The critical behavior close to the transition is that of an 
XY model in d-|-l dimensions. 

More interesting still is the behavior of as a function 
of the ns- The variation of is qualitatively different 
for the different parameter regimes: a) the weak coupling 
case for g < gmax, b) the intermediate one with gmax < 
9 < gcrit and c) the strong coupling limit for g > gent 
(see Fig. EJ. Going from the limit a to c we find that 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 

g/t 

FIG. 4: Phase diagram as a function of g/t describing the 
boundary between an insulating(I) and superconducting(SC) 
phase with long range phase coherence for the case of fully 
symmetric limit (A_b = 0,npi = npi = 2nB = 1). 




0.00 ' — — ' — — — — — — ' — — ' 

0.00 0.20 0.40 0.60 0.80 1.00 

FIG. 5: Phase boundary lines between the superconduct- 
ing(SC) and insulating(I) state at different values of the cou- 
pling and upon varying the average density of bosons. Solid, 
dotted and dashed line stand for values of the coupling con- 
stants for which g < gmax, gcrit > g > gmax and g > gcrit, 
respectively. 



with increasing ub the critical line decreases to zero, goes 
through a maximum starting at a finite T'^{nB = 1) and 
finally shows a SIT at a critical density. 



V. ROLE OF BERRY PHASE TERM: AN 
INTRINSIC MAGNUS FORCE ON VORTICES 

In the preceding section we have analyzed the basic 
physics resulting from the effective action by neglecting, 
(i) the influence of any feedback between the density and 
amplitude fluctuations (included in the field l[9i]) on the 
phase dynamics, and (ii) the contribution from the Berry 
phase term, responsible for the correct quantization of 
the pseudospin variable, which is given by the integral 
over all the possible paths of is{l — cos[0i])dr<j>i- In this 
section, we discuss the consequences of the presence of 
such a Berry phase term in the case where the phase 
action has a vortex solution. The existence of a vortex 
solution is assured for 2D systems where the phase cor- 
relations are described by a XY type dynamics. We shall 
show here that the Berry phase term will produce an in- 
trinsic Magnus Force on the vortex which is analogous to 
the Lorentz force for a charged particle, whose effective 
magnetic field depends on the spatial distribution of the 
6 field. Using the correspondence with the bosonic pseu- 
dospin variable, this implies a relation with the spatial 
density distribution of the bosons. The Magnus force has 
been widely discussed in the context of normal BCS type 
superconductors where it has been shown that it arises 
from the Berry phase caused by the adiabatic motion of 
a vortex along a closed loop coming back to its starting 
position^. The adiabatic vortex motion on a loop in the 
superconducting state turns out to be affected by an ef- 
fective magnetic field generated by the supercurrent aris- 
ing from the gradient of the phase which encircles such a 
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vortex. In the BFM scenario discussed here, we find that 
on top of the usual contribution due to the superfluid 
electrons, there is an intrinsic Berry phase term which 
will generate such a Magnus force and which has an in- 
tensity proportional to 2ns — 1, where ns is the bosonic 
density. This Magnus force, arising solely from the quan- 
tum nature of the pseudospin variable, is clearly indepen- 
dent on any superconducting state, and does not require 
a coherent superfluid current induced by the presence of 
the vortex itself. Moreover, since ub can be controlled 
externally, one has the possibility to tune the strength 
of the effective magnetic field acting on the vortex, and 
hence to alter its dynamical properties. This will result 
in a possible measurable effect on transport coefficients, 
such as the Hall coefficient, resistance, etc. 

To be more explicit and following the procedures used 
in different approaches treating with the Magnus force 
problem, let us assume that one has a vortex centered 
at the position R. Let us furthermore consider that 
we are in the continuum and at zero temperature so 
that the Berry phase term in the Lagrangian now reads 
Lb = J (Prcos[9{r,t)]dt4>{r,t). The contributions linear 
in dt(t>(r,t) and having time independent coefficients do 
not contribute to the dynamics. Moreover, since we con- 
sider that the phase 4> and the variable 9 are linked to a 
vortex solution centering at R, it is judicious to introduce 
a new relative variable r — R, which defines the intrinsic 
position dependence of the variables with respect to the 
vortex center. 

Starting from the above given expression for the La- 
grangian Lb and making a simple change in the time 
derivative we can rewrite it in the following form: 

Lb = J (frcos[e{r- R,t)]dt(j){r- R,t) (61) 
= — / (frcos[9{r - R, i)]V^(/)(r - R, t). (62) 

At this point, one recognizes that the term 
J <Prcos\9{r ~ R,t)]V r(t>{'r — R,t) plays the role of an 
effective vector potential As and ^ represents the ve- 
locity of the vortex. This Lagrangian is thus equivalent 
to that of a charged particle in presence of an effective 
magnetic field in the z-direction given by B = V x A. B 
hence creates a Magnus force Fm = ^ z)2Trcos[9o] 

which is transverse and proportional to the vortex ve- 
locity and whose strength is linked to the magnitude of 
this effective magnetic field = —2Trcos[9o]. Its magni- 
tude is given by the asymptotic value of the background 
bosonic density at large distance cos[0o] for r — s- cx). 

In case of a 2D superconductor described by the BFM 
scenario, the usual effective magnetic field —2'Kps aris- 
ing from the supcrcurrcnt of the electrons circling the 
vortex has thus to be supplemented by the above men- 
tioned contribution B (intrinsic to such a scenario) and 
permits in principle to modulate the total Magnus force 
upon changing the density of the bosons. Clearly, a full 
description of such an eventuality, which is beyond the 
scope of the present analysis, has to be studied in more 



detail and has to include the derivation of the effective 
action in presence of magnetic field such as to determine 
the superfluid density of the electrons ps- 

At this stage we simply want to point out that the sign 
of the magnetic field arising from the Berry phase term 
in the action of the BFM changes if one tunes the bosonic 
density between zero and unity. In terms of the variable 
cos[^] this implies a variation in the range [—1,1]. In 
other words, the sign of the intrinsic Berry phase induced 
Magnus force will change at ub = 1/2 when going from 
the limit of small density of bosons to high density. This 
extra contribution to the Magnus force has to be added 
of course to the conventional one known for standard 
superconductors. 



VI. COMPARISON WITH NEGATIVE-U AND 
BOSE-HUBBARD MODELS 

A frequently asked question concerns the qualitative 
differences which exist between the BFM and similar sce- 
narios such as the negative-U Hubbard model and pure 
bosonic systems such as the Bose-Hubbard model and 
Josephson junction arrays. 

Let us start with a comparison of the BFM and the 
negative-U Hubbard model which has been studied in a 
great variety of different approaches and discussed espe- 
cially in connection with the BCS-BEC crossover. The 
main issue of such a comparison is the occurrence or not 
of a quantum critical point in the negative-U Hubbard 
scenario at a finite value of the coupling U which intro- 
duces the pairing among the electrons. As we have shown 
in this paper, in the BFM there exists a critical value for 
the exchange coupling g/t responsible of the pairing in 
the fermionic subsystem above which the system under- 
goes a SIT. We have seen that within a path integral 
formalism, one can extract an effective action describing 
such a transition as being primarily due to the interplay 
between the local quantum pair exchange (between a bo- 
son and a fermion pair) and the non-local inter-site pair 
hopping leading to a phase dynamics originating from the 
itinerancy of the fermions. Due to the non-trivial depen- 
dence of the strength of these two competing mechanisms 
on the microscopic parameters of the BFM, there occurs a 
SIT for a finite density as well as finite exchange coupling 
g/t. For the negative-U Hubbard scenario this is not the 
case and the SIT does not occur at any finite value of the 
ratio U/f^^. This model merely describes a continuous 
cross-over between a BCS superconductor and a BEC of 
tightly bound electron pairs as U is increased from to 
oo. 

In order to better understand this difference between 
the BFM and the negative-U Hubbard model, let us con- 
sider two scenarios in equivalent situations, i.e., the half 
filled band case for the negative-U Hubbard model and 
the fully symmetric case for the BFM (with the bosonic 
level lying in the middle of the fermionic band such that 
both, the fermionic band as well as the bosonic level are 
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half occupied). We then address the question how, in the 
strong coupUng hmit, the pair hopping is generated out 
of the basic configurations of states and processes which 
contribute in this regime. 

For the case of the negative-U Hubbard model with 
fermions interacting via a local attractive potential, the 
ground state wave function in the limit of J7 — s- cx) and 
at half-filling is highly degenerate and is composed of all 
the possible configurations comprising equal distribution 
of zero and doubly occupied site. In this limit, one can 
perform a mapping of this model on the hard-core boson 
model described by the Hamiltonian: 

i/£/ = -/i^n,. (63) 

bf {bi) and n.i — b^bi stand for the creation (destruc- 
tion) operators of hard-core bosonic particles (tightly 
bound electron pairs) and for their density operator, re- 
spectively. Due to the presence of a coupling, which is 
isotropic both in the boson hopping and in the charge 
interacting channel, one has a superconducting state for 
any finite value of the ratio U/t, and possibly a super- 
solid phase due to the symmetry of the charging interac- 
tions characterized by a coexistence of diagonal and off- 
diagonal long range order. This implies that the quantum 
critical point is strictly pushed to U —foo. 

In the BFM the phase space in the large g/t limit is 
completely different. In the fully symmetric case of this 
model (corresponding to a total density equal twice the 
number of sites {ritot — 2)), the ground state configura- 
tion is given by a wave function where all the bosons are 
strongly coupled with pairs of fermions, thus resulting in 
a product state of local bonding states: 

This wave function can be viewed as a ferromagnetic 
Ising-type state in the sense that it is made up of a ferro- 
type bonding order and is not degenerate with respect 
to all the other possible local configurations. Moreover, 
by construction it contains bonding-bonding correlations 
on a long range (bonding solid), but has zero phase cor- 
relation length. The latter can be shown by evaluating 
the static correlation function for the bosons or, equiv- 
alently, the fermion pairs which gives {tpo\pt Pjli^o) = 
and (ipolT^Tjlipo) = for any distance |i — j|. 

Next, let us consider the low energy configurations 
which are mixed into such a, g/t = oo ground state when 
the kinetic energy operator Ht = Y^i^j, cr^iji^ta'^ja + 
H.c ) is switched on. Applying Ht to li/'o); there will 
occur states which are separated by an energy gap with 
respect to |i/'o)- The low energy states to be consid- 
ered as the relevant quantum mixed configurations are 
of non-bonding nature, such as Cf^\Q) = c^|0) and 
•^ilrlO) — Pi^tA^)- order to construct those states, let 
us to begin with consider a local excitation on two ad- 
jacent sites i,j given by Sf^C'^_^\Q) in the background 



of bonding states. This configuration can be seen as a 
ferro-type order interrupted by two domain walls of non- 
bonding/bonding nature. Now let us consider the dy- 
namics of such objects and under which conditions they 
can be rendered itinerant in a way which leaves the num- 
ber of bonding configurations unchanged. The problem 
is thus analogous to that of domain wall dynamics in an 
Ising type systems. It so happens (a detailed discussion 
of this is beyond the scope of the present study and will 
be given at a later stage elsewhere) , that the local degen- 
eracy of the non-bonding states will be removed by the 
action of the kinetic term. In this way it induces a global 
lowering of their excitation energy due to the dispersion 
of each domain wall which is of the form 



Edw (^) = 2 ["'^ + ^ v/(efc - 25)' + 4<2] (65) 

and where e(fc) = —2ztcosk. Only after the condensation 
of the domain wall like excitations in the presence of the 
bonding state background one can meet the conditions 
for setting up long range boson and fermion-pair phase 
correlations . 

This simple sketch of the nature of the excitations in 
the BFM implies that the onset of phase coherence can- 
not be activated for an arbitrarily small hopping am- 
plitude since one has to overcome the energy gap be- 
tween the ground state and the manifold of non-bonding 
states, which is of the order of ~ 2g. As we can see from 
the expression for the dispersion of the domain walls, 
this is achieved when i?£)n/(fc) becomes zero for fc = 0, 
i.e., when g/t — |;(2 + \J2^ Az'^). For z=2 this gives 
g/t ^ 2.0, which determines when those defects-like ex- 
citations become gapless and thus induce a proliferation 
within the background of the bonding states. It is worth 
pointing out that this value for the critical exchange cou- 
pling matches rather well the value obtained from our 
functional integral approach discussed in this paper (see 
Fig 3). ^ 

The picture which hence emerges for the BFM is very 
similar to that of a XY model in a transverse field. There, 
the hopping term is responsible of the XY dynamics, 
while the local boson-fcrmion pair exchange provides the 
role of the effective transverse field. This itself is already 
a strong indication that this model is of different nature 
to that of the negative-U Hubbard scenario, for which 
we know that it is, on the contrary, akin to an isotropic 
Heisenberg model. 

A further essential difference between the BFM and the 
negative-U Hubbard model can be seen from their respec- 
tive path-integral formulations. When the path integral 
representation is considered for this model^, the first 
step is to make use of a Hubbard-Stratonovich transfor- 
mation to rewrite the quartic interaction term in a bilin- 
ear form, where the fermionic operators are now coupled 
to random auxiliary fields. Limiting oneself to purely 
superconducting order, the usual procedure for manipu- 
lating such a bilinear action is to separate the complex 
auxiliary field into its modulus and a pure phase part 
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before performing an expansion around the saddle point 
solution for the amplitude of this auxiliary field. This 
way of proceeding allows then to extract an effective ac- 
tion for the slow phase dynamics. This is distinctly dif- 
ferent from the functional integral representation for the 
BFM which we have presented above and where from 
the very beginning and throughout such a procedure the 
fermionic pair fields are coupled to physically real bosonic 
modes which have their own proper dynamics. As we 
have seen in section III, the bosonic part of the action 
is treated within a coherent pseudo-spin representation 
where the dynamics of the pseudospin is parameterized 
by the time dependence of the spherical variables. The 
role of the dynamics within the spherical representation 
of the bosonic field and of the Berry phase term (which 
is a consequence of the quantum interference in the local 
pseudospin space) is a distinct feature of this BFM and 
presents specific differences with respect to the negative- 
U model. Such differences imply in particular that in the 
BFM case one has feedback effects between the ampli- 
tude and the phase fluctuations which are totally absent 
in an equivalent description of the negative-U Hubbard 
model. Such features are important because of intrinsic 
dissipation effects that eventually can change the nature 
of the transition and possibly are relevant for the emer- 
gence of a "bosonic metal" in proximity of the SIT, a 
feature which seems to be outside the framework of the 
negative-U Hubbard model. 

Let us conclude this comparison of the BFM with sim- 
ilar models with a brief discussion on the Bose-Hubbard 
model and Josephson junction arrays. The main aspect 
which emerges as a common denominator in all those 
models, at least as far as the phenomenon of the SIT 
is concerned, is that the mechanism which is responsi- 
ble for the degrading of the phase coherence is analogous 
and originates from the competition between the phase 
and charge degrees of freedom. In the Bosc-Hubbard 
model this manifests itself as a competition between the 
boson hopping and their charge repulsion, while in the 
Josephson junction array scenario it appears as an in- 
terplay between the Josephson tunnelling amplitude and 
the charging energy. 

In spite of this, at first sight, similarity between the 
BFM and those scenarios, we would like to stress that in 
a system where the dynamics is described by the BFM 
scenario, the interplay between the phase and amplitude 
fluctuations is intrinsically related to the coupling be- 
tween the fermionic and bosonic degrees of freedom. This 
introduces amplitudes of the different processes at work 
in form of a non-trivial dependence on the microscopic 
parameters of the starting Hamiltonian. Furthermore, 
as we have seen in the above discussion of the effective 
action of the BFM, this model contains features which 
go beyond those which characterize the pure Joscphson- 
type dynamics and which arise from the peculiar feedback 
between the fluctuations of the bosonic density (or am- 
plitude pair field) and fluctuations of the phase. Last not 
least, the appearance of an intrinsic Berry phase-term in 



the BFM scenario together with the cflect of dissipation 
due to the fermionic dynamics, can give rise to an uncon- 
ventional phenomenology when topological phases play a 
role, and especially in the presence of vortices. 



VII. CONCLUSION 

In this paper we examined the nature of a superconduc- 
tor - insulator transition in a system of localized bosons 
and itinerant fcrmions coupled together via a pair ex- 
change term. An effective action was derived from such 
a microscopic model, which, after integrating out the 
fermionic fields, could be phrased in terms of amplitude 
and phase fluctuations of the bosons. In order to make 
the presentation more familiar we discussed the action in 
a terminology frccjucntly employed in connection with the 
study of Josephson junction arrays. We stress that our 
system does not necessarily imply any charged fermions 
and bosons. 

Considering the phase-only part of the effective action 
it is fully equivalent to the quantum phase model for the 
Josephson junction arrays, discussed in terms of: (i) a 
Josephson coupling term, (ii) a charging or capacitance 
term and (iii) an offset charge term. Equivalent to that 
in our scenario is: (i) a boson hopping term, (ii) a term 
which takes into account the reduction in hopping am- 
plitude due to a fluctuating local boson density arising 
from the intrinsic on-site exchange coupling between the 
bosons and the fermions and (iii) a chemical potential 
term controlling the bosonic concentration. New in the 
present study within already the lowest (phase-only) ap- 
proach to the boson-fermion system is the intricacy of 
the dependence of the effective Josephson coupling, the 
capacitance term and the off-set charge term on the pa- 
rameters of the initial Hamiltonian, i.e., the exchange 
coupling and the density of bosons (or else the total par- 
ticle density). It turns out that within the lowest (phase- 
only) contribute a superconductor to insulator transition 
can not only be triggered by a change in the exchange 
coupling but also by a variation of the boson density. 
The latter presents evidently a certain interest from the 
experimental point of view and can possibly be tested in 
such transition occurring in optical lattices for ultracold 
fermi gases with Feshbach resonance pairing. 

Apart from the phase-only part of the effective ac- 
tion we established the existence of an intrinsic Berry 
phase term, which arises from the hard core nature of 
the bosons and gives rise to an additional Magnus force 
when the system is such that topological ground state 
configurations like vortices are stabilized. This again is 
of potential interest for experiments since in principle one 
can change the sign of the Magnus force upon changing 
the density of bosons and thus deviating the motion of a 
vortex from one direction into the opposite one. These 
very preliminary results will be dealt with in greater de- 
tail in some future studies. 

A further property of the Berry phase term is that it 
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gives rise to a bilinear coupling term between the phase 
and amplitude fluctuations, and is hence much more di- 
rect and relevant than similar terms in, for instance, sce- 
narios based on the negative-U Hubbard model where 
they occur only at a much higher order in a correspond- 
ing loop expansion of the trace in the effective action. 
This again merits to be investigated in some detail with 
the aim to study the dissipation introduced by such am- 
plitude phase coupling and its effect on the nature of 
the transition, in view of exploring the possibility of an 
intermediary bosonic metallic ground state. 

Finally, we have dealt with a frequently asked question 
concerning the differences between the negative-U Hub- 
bard scenario - mainly studied in connection with the 
BCS-BEC cross-over and the presently studied boson- 
fermion model. Far from being able to give a complete 
account for the major differences we found mainly two 
aspects which distinguish the physics of these models on 
a qualitative and robust level. The one is that in the 
negative-U Hubbard model a superconductor - insulator 



transition can not take place at any finite coupling U, nor 
can such a transition be triggered by the change in par- 
ticle concentration. The second point is that the ground 
states in the strong coupling limit of the two models are 
quite different: a highly degenerate ground state for the 
negative-U Hubbard model with excitations being con- 
trolled by an isotropic Heisenberg model when the hop- 
ping term is switched on. Contrary to that, for the boson 
fcrmion model the ground state is non-degenerate (cor- 
responding to a ferro pseudo-magnetic Ising type system 
of singlets formed by bosons and fermion pairs), which, 
after switching on the fermion hopping term, gives rise 
to propagating domain like structures. The topological 
structures appearing in the ground state might have mea- 
surable consequences in the transport properties near the 
superconductor-insulator transition. 

This and the other preliminary studies mentioned 
above will require a detailed analysis and will be dis- 
cussed in future. 
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